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Abstract 

A bacterial colony may develop a small number of cells genetically 
identical to, but phenotypically different from other normally growing 
bacteria. These so-called persister cells keep themselves in a dormant 
state and thus are insensitive to antibiotic treatment, resulting in se- 
rious problems of drug resistance. In this paper, we proposed a novel 
strategy to "kill" persister cells by triggering them to switch, in a fast 
and synchronized way, into normally growing cells that are susceptible to 
antibiotics. The strategy is based on resonant activation (RA), a well- 
studied phenomenon in physics where the internal noise of a system can 
constructively facilitate fast and synchronized barrier crossings. Through 
stochastic Gilliespie simulation with a generic toggle switch model, we 
demonstrated that RA exists in the phenotypic switching of a single bac- 
terium. Further, by coupling single cell level and population level sim- 
ulations, we showed that with RA, one can greatly reduce the time and 
total amount of antibiotics needed to sterilize a bacterial population. We 
suggest that resonant activation is a general phenomenon in phenotypic 
transition, and can find other applications such as cancer therapy. 
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1 Introduction 



Noise has often been viewed as a nuisance for many years in biofogy. Robustness 
of large biological systems requires noise from both intracelluar and intercelluar 
sources being canceled or filtered in one way or another. Yet growing evidence 
indicates that noise actually plays fundamental roles in many biological pro- 
cesses, as in cell fate decision and in mutation and evolution [T|. In this paper, 
we suggest that the functional role of noise can go even beyond; the internal 
noise in bacterial gene expression can be utilized to counteract antibiotic resis- 
tance, by inducing resonant activation that can facilitate a fast and synchronized 
phenotypic switching in bacteria population. 

Antibiotic resistance is a severe and growing problem in clinical practice. It 
refers to a phenomenon that certain phenotypes of microorganisms, e.g. bac- 
teria, are able to withstand (and requires prolonged) antibiotic treatment. It 
may be acquired from horizontal gene transfer and mutations in the pathogenic 
chromosome [H 12111] , or from the existence of phenotypic heterogeneity within 
bacteria population [5J [51 [7] • The latter links antibiotic resistance to a special 
bacterial phenotype called persister cell, a non-growing (or slowly growing) and 
non-inherited cell phenotype whose number only accounts for a small fraction 
of total population. 

Persister cells are genetically homogeneous to normally growing cells. They 
are first discovered by Bigger "8^ and then found in higher-percentage in biofilms 
that are known to be responsible for a majority of recalcitrant infections such 
as tuberculosis [HI HH]- Several experimental and theoretical works suggested 
their formation during mid-exponential phase, as well as their function as an 
"insurance" to maximize the overall survival probability of bacterial population 
in changing environments [TTl [T^ [T3] . The magic is that persister cells have 
negligible growth rate and non-negligible phenotypic switching rate. The negli- 
gible growth rate helps persister cells dodging antibiotic attack which depends 
mostly on active cell wall growth. The non-negligible phenotypic switching rate, 
on the other hand, ensures finite probability of stochastic switching from per- 
sister cells to normally growing cells taking place at a period that the stress 
(e.g., antibiotics) is removed. Then those newly formed normal cells serve as 
the "seeds" for reestablishing the population. Therefore, Kussell et al. proposed 
the phenotypic switching rate can be seen as a result of evolutionary adaptation 
of bacteria to their real fluctuating environment |12| . 

The ubiquity of persister cells makes bacteria population hard to sterilize. 
The time-series of the survival fraction under antibiotic treatment obeys a two- 
phase exponential decay, with the majority of cells being killed at a fast rate at 
the beginning while the rest being killed much slowly afterwards. hipA7, a mu- 
tant strain of E.coli that contains higher percentage (10~^ ^ 10~^) of persister 
cells, has been reported survival in a fraction about 10^^ even after continuous 
ampicillin treatment for 50 hrs [5]. Therefore, It is of concern how to efficiently 
sterilize bacteria populations, especially for strains with more fraction of per- 
sister cells. In 2008, Gefen et al. observed that persister cells of hipA7 assume 
normal growth during the first 1.5 hour on exiting the stationary phase |14| . 



While applying ampicillin within that time window, the number of persister 
cells can be lowered by 1 order of magnitude. However, those still alive may 
adopt dormancy after that time window, and can convert to normally growing 
phenotype and re-grow to a new population under appropriate conditions. The 
essential problem here is that well-established persister cells are insensitive to 
antibiotics. They have to be converted into normally growing cells to get steril- 
ized by drugs. However this transition is stochastic and may take a long time. 
For hipA7, the transition rate is 0.07/ir~^, which gives an exponential waiting 
time distribution with the average waiting time of conversion 14 hr. For some 
species, e.g. the E coli mutant hipQ, the rate can be much smaller. Therefore, 
bacteria sterilization requires continuous antibiotics application at least to cover 
this broad range of time, which may be impractical and/or detrimental to the 
host. 

Several toxin-antitoxin (TA) modules in bacterial chromosome have been 
experimentally identified to regulate bacterial phenotypic transitions, though 
detailed mechanism has not been clearly understood(71 [T51 [TB]. Irrespective of 
the details, the basic mechanism is simple: a mutual inhibition exists between 
antitoxin's and cognate toxin's expression, which determines whether a single 
bacterium assumes normally growing phenotype (when antitoxin dominates) 
or persister phenotype (when toxin dominates). Therefore bistability is the 
major dynamical property of single cell's phenotypic transition [T^. In this 
theoretical investigation we used the well studied generic toggle switch to model 
this bistable system. Figure 1 (a). 

Making analogy between biological networks and other familiar physical sys- 
tems has lead to several insightful studies [ITl [HI [19] . In this work we focused 
on resonant activation (RA) , a well studied phenomenon for thermally activated 
barrier-crossing systems [2TJ [35] . If the barrier is under time varying periodic 
perturbation, "cooperative interplay between the barrier modulation process 
and thermal noise assisting barrier crossing events can cause an enhancement of 
the reaction kinetics" [20]. The mean first passage time (MFPT), which is the 
average time the system waits for the first successful barrier crossing, reaches its 
minimum (by several folds or even orders of magnitude compared with that of 
the unperturbed system) at the resonance frequency of the perturbation. Under 
resonance frequency Uc with period T^^^ = 27r/wc, the system prefers to make a 
transition when the barrier height reaches its minimum. Consequently, the FPT 
distribution displaces a series of peaks at odd multiples of 11^^/2, instead of a 
continuous exponential distribution. RA is related but different from another 
well studied phenomenon, stochastic resonance [201 [23] ■ 

Complimentary to current efforts of searching for more efficient antibiotics, 
in this work we propose to utilize the phenomenon of resonant activation to 
help fighting drug resistance. Noticing the similarity between thermally acti- 
vated barrier crossing and cell phenotypic transition, we will first use a genetic 
toggle switch model to demonstrate that resonant activation exists for biological 
network dynamics. Then we will examine how one can shorten the time and the 
amount of antibiotics needed to extinct a bacteria population using resonant 
activation. Our strategy utilizes the two characters of RA: accelerated kinetics. 



and multi-peaked FPT distributions. 



2 Model and Method 
2.1 Single cell level 

Currently the exact regulation mechanism for the persister-normal cell transition 
is not clear, and different hypotheses have been raised to give possible answers[7]. 
For our purpose of illustrating the idea, we will use a generic toggle switch to 
represent the mutual inhibition within T-A module. As shown in Figure 1 
(a), the network contains two genes mutually inhibiting each other through 
their dimerized protein products. The two types of protein dimers compete 
for the promotor binding site. Ten chemical reactions shown below control the 
dynamics of the generic toggle switch shown in Figure 1 (a). 

0%0 + A,0%0 + B 

k, Jcr 

A + A^A2,B + B^B2 

+ A2^ OA2, O + B2 ^ OB2 

OA2 H OA2 + A, OB2 OB2 + B 

Table 1 gives the corresponding rate constants in reduced unit. With these 
parameters and mass-action type dynamics, the system gives two stable states 
corresponding to the persister (the difference between the total number of free 
and bound protein B and the total number of free and bound protein A, denoted 
as (5, is less than 0), and normally growing {5 > 0) phenotypes, respectively (see 
Figure 1 (b) ) . Stochastic fluctuations of the reactions drive the system transit 
between the two states. These parameters are modified from the model of Allen 
et al. [21] , so the model gives hip A 1 p2n (persister cell to normally growing cell) 
and n2p (normally growing cell to persister cell) switching rates 0.07/ft.r and 
0.008//ir (with time unit 1 unit = 130 hr~^), respectively, as used by others [3 
\n\ . The abstract toggle switch model also places the present work in a broader 
context. Many microorganisms, including viruses and bacteria, coexist in a 
dormant and an active phenotype |25j . The toggle switch is a frequently occured 
generic regulation mechanism and a good model for phenotypic transitions in 
bacteriaf^ITT]. 

To be specific and for practical considerations of the computational feasi- 
bility, we chose model parameters in most simulations in this work to mimic 
the dynamics of hipAl . However, we want to emphasize that the proposed ap- 
proach below works best for the following situation. First, we assume that one 
can regulate some of the rate constants through an external oscillating perturb- 
ing signal. Consequently some rates (the degradation rate of A for the results 
reported in the main text) are oscillating with time. The strength of the per- 
turbing signal should be restricted due to the consideration of toxicity to the 



host. Second, we focus on the case that the transition from the persister to the 
normahy growing phenotype is very slow, so that the sterihzation of the total 
population requires a long term antibiotic treatment, which may bring severe 
side-effect like liver damage. 

2.2 Population level 

At the population level, a cell is subject to an environment alternating between 
growing and antibiotics stress conditions. A normally growing cell has faster 
net proliferation rate than a persister cell does under growth condition, but also 
a larger death rate under antibiotic stress (see also |I2j). Three types of cellular 
events can take place for each cell: cell division-the cell including the molecular 
state of the goggle switch are cloned into two identical copies, cell death, and 
phenotypic transition. Table 2 gives all the related rate constants based on 
experimental observations [5l [12] . For a given cell with phenotype n (normally 
growing) or p (persister), the reactions for stochastic simulations additional to 
the 10 toggle switch reactions are, 

n ^ 2n, n'^ (p, p 2p, p (f). 

Unlike the work of Kussell et al, in our model we did not simulate the cell 
phenotypic transitions directly. Instead we propagate the 10 toggle switch reac- 
tions for each cell, which determine the phenotype of the cell (see above). Each 
population level simulation initiates from a stationary-phase colony including 
10"* normal cells and 10^ persister cells. The population is then put into fresh 
medium with/without antibiotics, as determined by each different strategy. The 
phenotype of a single cell is determined by its own toggle switch dynamics based 
on the value of 6. We use r-leap Gillespie algorithm to propagate the 10 chem- 
ical reactions of each toggle switch and cell birth/death simultaneously [28] . 
For the rates with periodic time-dependence, we approximate them as constant 
within one Gillespie step, which is much smaller than the rate varying period. 
For simplicity, we do not consider quorum sensing, thus each cell behaves inde- 
pendently and does not communicate with others except for competing resources 
as discussed below. To prevent the population from overgrowth in our simula- 
tions, we rescaled the growth rate g as a decreasing function of the number of 
normally growing cells n, gn{t + t) — .gfj/[l + a ■ n{t)] where r is the time step 
in r-leap Gillespie algorithm, a controls the scaling strength (here a = 0.001), 
and g^i^ represents the original growth rate of normal cells without rescaling. We 
didn't rescale the persister cells' growth rate since the value before rescaling is 
already negligible. A physical justification of the rescaling is that the accessible 
nutrients and volume of a cell colony are usually limited against unrestricted 
massive replications. 



3 Results 



3.1 Resonant activation exists in cellular phenotypic switch 

For the model we examined, the switching rate from n2p is much smaher than 
that of p2n at the exponential phase. This dynamics mimics that of hip A 7 [5] , 
and resembles to barrier crossing rates in an asymmetric double well potential, 
with one deeper well representing normally growing phenotype, and the other 
well representing persister phenotype. Let's define ip2n as the first passage time 
for p2n switching. The distribution of tp2n, apart from an initial transient time 
period, follows an exponential form 

P{tp2n) OC exp{-tp2n/tK) 

(see Figure 2 (a)), where tK is the inverse of the Kramers rate [23l [29l [30] . The 
prolonged distribution contributes to bacterial persistence. 

Next, we perturbed the protein A's degradation rate with sine- formed signal 
(see Figure 1), = fcg [1+7 sin(wt-|-0)], with < 7 < 1. For each simulation the 
phase 9 is randomly drawn from a uniform distribution betwen and 2tt. This 
is because in the population level simulations below, the relative phase between 
the birth time of a persister cell and the added signal can be seen as a random 
variable. Figure 2 (b) shows that the mean first passage time (MFPT) as a 
function of w, averaged over 5000 independent simulations, shows a minimum 
around Wc = 1.6 /ir~^, although the curve is rather fiat over a range of ui values. 
Existence of the minimal MFPT signatures resonant activation [20l [22] , despite 
the current system is described by discrete dynamics. 

Figure 2 also shows P(ip2n) under the perturbing signal with different fre- 
quency. At a high w, the system cannot respond fast enough, and the perturba- 
tion is equivalent to an averaged constant one. The distribution is exponential, 
Figure 2 (c). Under resonance frequency Wc, however, tp2n distribution changes 
into several separated spikes. Figure 2 (d). Note that the peaks of spikes over- 
lap with the peaks of the periodic signal. Therefore transition takes place more 
frequently when the signal reaches the peak value, thus the degradation rate of 
protein A is the fastest, and the transition from A dominant to B dominant is 
the easiest [23]. If the system misses one peak of the signal for a p2n transi- 
tion, it prefers waiting for the next peak. We observed this type of localized 
distribution over a broad range of w values corresponding to the fiat bottom 
region of the MFPT-w curve (see Figure 2 (b) and Figure 6 in the Supplemen- 
tary Information). Further increase of to leads to gradual merge of the spikes, 
and eventually reduce to a single exponential distribution. These observations 
are consistent with studies with barrier crossings in a double well potential [20] , 
further supporting the existence of RA in the system. 

For all the results reported here, we added the perturbing signal on protein 
^'s degradation for illustrative purpose (see Figure 1 (a) ). In real experiments 
and applications, one can also perturb other reactions {e.g. protein synthesis as 
well as degradation), depending on the actual practical feasibility. For example, 
one possible implementation of this perturbation may be through varying the 



activity of protease through specific regulating molecules. We further presented 
results with perturbations on A's synthesis rate, and B's synthesis and degra- 
dation rates in the Supplementary Information, Figure 7. In all these cases we 
observed resonant activation with the same resonant frequency, but with vary- 
ing fold of change of the MFPT. We provided a theoretical explanation there. 
It depends on the system to identify the reactions most sensitive to the pertur- 
bations. Figure 8 in the Supplementary Information also showed that existence 
of RA is independent of the detailed form of the periodic perturbation signal. 

3.2 Resonant activation accelerates bacteria colony ster- 
ilization 

The stochastic simulations on single cell dynamics discussed above show that 
resonant activation can facilitate fast and synchronized p2n switches. Next we 
coupled the single cell level dynamics with population level proliferation/death 
under changing environment. 

First we define a killing strategy K to be the combination of a perturbing sig- 
nal S and an antibiotic environment E, characterized by their adding frequency 
(wi and UJ2), strength, and duration, respectively. 

For a given population initialized with 10^ normally growing cells and 10^ 
persister cells, our simulations allow it to evolve until no bacterium exists, or a 
maximum time reaches. Because the stochastic nature of the dynamics, for each 
population the sterilization time T^m is random. To quantitatively compare 
different strategies, for each strategy we performed independent simulations with 
1000 populations, and recorded Tkui of each population sample. Figure 9 in the 
Supplementary Information gives the killing time distribution of 1000 population 
samples under strategies Kl and K3. To compare the different strategies, we 
used the time needed to sterilize 90% of the 1000 population samples, Q, as a 
criterion. Practically this is a more relevant quantity than the average killing 
time, although we reached similar conclusion below with the latter. 

The stochastic simulation results are summarized in Figure 3 and Table 1. 
Figure 4 give several typical trajectories. Without the perturbing signal, both 
the strategies with periodic (Kl) and continuous (K2) antibiotics treatment re- 
quire long time, since it takes long time to eradicate the persister cells. On the 
other hand, with the perturbation at the resonance frequency (K3 an K4), the 
sterilization time is greatly reduced. Because under K3, most of the p2n tran- 
sitions take place within the period of applying antibiotics, it is difficult for the 
bacteria population to restore cither the normally growing or the persister sub- 
populations (see Figure 10 in the Supplementary Information). Consequently 
the sterilization time for K3 and K4 arc about the same. To further prove this, 
we compute over 500 independent samples the ratio {Rp2n) between the number 
of p2n transitions that happen during growth environment and the total number 
of p2n transitions. Larger i?p2n corresponds to the inefficiency of the periodic 
antibiotic strategy, because the population may be easier to get re-estabished 
during growth period. Under K3, R^2n — 4-3%, while under Kl, R^^n ^ 46.8%. 
We observed similar efficient bacterial eradication with signal frequencies away 



from u)c but still laying near the flat bottom of the MFPT w curve in Figure 1 
(b). However, further change of the frequency (K5 and K6) shows less improve- 
ment over that of Kl or K2. Figure 3 (b) gives the total amount of antibiotics 
used for each case. Compared to K2, K3 requires less than half of the time, 
with ~ 25% of the total amount of antibiotics. In the Supplementary Informa- 
tion, we also examined how the killing time depends on the signal duration and 
strength, Figure Sll and Figure 12. 

To further examine the strength and limitation of our proposal, we exam- 
ined a more difficult case of bacterial persistence. The model is similar to what 
discussed above, except that the p2n transition rate is reduced by a factor of 
5. Correspondingly, we extend each period of alternative antibiotics and freely 
growing environment. This leaves more time for the system to recover under 
the latter environment. Figure 5 gives several typical trajectories of population 
level simulations. For strategy Kl, we had difficulty to observe population ex- 
tinction even after 10'' hours. While almost all the normally growing cells are 
killed, and the persister subpopulation size is reduced under an antibiotics en- 
vironment, both subpopulations are restored to their original levels under next 
freely growing period. For strategy K2, the 90% quantile of sterilization time 
for 1000 populations is 362 hrs. In comparison, the 90% quantile of sterilization 
time for K3 and K4 are 358 hrs and 168 hrs, respectively. Though under K3, 
the ratio of p2n transitions within growth period Rp2n — 4.3%, in this case a 
single p2n transition under the freely growing environment is sufficient to re- 
store both the two subpopulations. That explains why K3 requires almost the 
same sterilization time with K2. However, K3 uses about half of the amount 
of antibiotics needed under K2. In addition, since antibiotics are applied pe- 
riodically, K3 might be better than K4 considering it leaves time for the host 
to recover from possible side effects of the antibiotics treatment. Therefore, in 
this case, one has to make a compromise between short sterilization time and 
side-effect reduction. 

4 Discussions and concluding remarks 

In the past decade, bacterial persistence became a spotlight in mircobiologi- 
cal arena. The leading actors, called persister cells, are some special dormant 
cells which account for only a small fraction of the total population. How- 
ever, they are insensitive to antimicrobial therapy, and are able to switch back 
into normally growing phenotype to initiate population regrowth. The pheno- 
typic switching property of bacteria has been experimentally identified to be 
regulated by several toxin-antitoxin (TA) modules within bacterium chromo- 
some, such as hipBA, relBE and chpA. For example, Keren et al. reported that 
/ijp_B^-knocked-out E. coli biofilm produced 150-folder-fewer persister cells un- 
der mitomycin treatment |31| . Continuous efforts have been made to reveal the 
detailed molecular regulation mechanism [151 116| . 

Now, it is believed that persister cells may evolve into maximizing the overall 
survival possibility of bacteria population in real fluctuation environment. Many 



lab experiments have proved that the existence of persister cells in biofilms is 
responsible for many recalcitrant diseases, such as human tuberculosis, an infe- 
cious disease caused by Mycobacterium tuberculosis biofilms. The typical antibi- 
otic treatment of this disease is as long as 6-9 months [3]. Therefore, besides the 
inefficiency of the therapy, the side-effect from such a long-term use of antibi- 
otics is of serious concern. Active research is undertaken to fight against bacteria 
persistence through accelerating the phenotype transition rate with chemical or 
physical method [71 [TU]. In this work we assume such a mechanism exists, and 
focused on the optimal strategy to combine it with antibiotics treatment. 

Stochastic resonance has been discovered in many biological systems [5^ . 
Similar to Schmitt et al., we also observed stochastic resonance in the present 
system. Resonant activation, on the other hand, is seldom discussed in the 
biology context. The phenomenon of RA is related to fluctuation resonance 
previously reported j33l |34]. For a stochastic system with one steady state 
driven by an oscillating perturbation, there may exist a resonant frequency of 
the perturbation so that the system shows largest fluctuation amplitude. This is 
called fluctuation resonance, which resembles the peak of the energy absorption 
spectrum of a system. For a system with multiple steady states, larger fluc- 
tuations lead to faster transition rate to a new steady state. This is resonant 
activation. Here we demonstrated the existence of RA in non-thermal systems 
governed by discrete chemical reaction dynamics. RA has two unique properties: 
reduced mean first passage time, and localized transition in time. We propose 
to utilize the two properties of RA to help on eliminating persistent bacteria. As 
previously mentioned, the essential reason for bacterial persistence is the large 
time-scale of tp2n distribution. This requires continuous antibiotic treatment to 
cover most of the time period to prevent bacteria population re-establishment. 
With RA, the distribution is narrowed through reduction of the MFPT. This 
is our main argument for using RA against bacteria persistence. Furthermore 
for some cases the localized spike shaped transition time distribution may al- 
low dividing the antibiotics treatment into sessions without serious problem of 
bacteria population restoration. This is of special advantage by minimizing side 
effects of antibiotics treatment to patients if the treatment has to be long. On 
the other hand, our simulations show that even if one can accelerate the bac- 
teria phenotype switching rate, improper procedure (strategy K5) leads to no 
improvement in fighting bacteria persistence. 

In this work, we presented the general idea of modulating cellular pheno- 
type switching through resonant activation. It should be viewed as illustrative. 
More detailed modeling and experimental studies are necessary to examine the 
feasibility and the optimal strategy for each specific system. Detailed molecu- 
lar mechanism of the toxin- antitoxin module, its interaction with related signal 
transduction and metabolic pathways should be carefully considered. The de- 
tailed model will provide information on the strategies of adding the perturbing 
signal. In a bacteria colony, several persistent phenotypes may coexist. In this 
case a more efficient strategy may be to apply antibiotics continuously for a 
period, then switch to the strategy we propose here for the most persistent phe- 
notype. The broad range of resonance frequency of one phenotype (see Figure 



2 (b)) may also allows one to choose an overlapping frequency for all the phe- 
notypes. The perturbing signal is not limited to chemicals, but any external 
environmental change that can affect the phenotype switching dynamics. 

We want to point out that the actual performance of each strategy depends 
on the property of the system, especially the unperturbed p2n transition time 
and the system noise level. Here we only focused on illustrating the basic idea, 
and the choice of the system is partly restricted by computational considerations. 
Orders of magnitude reduction of the MFPT with RA have been reported for 
some physical systems [3S] . One may expect similar result for some phenotypic 
transitions. 

While here we focused on bacterial persistence, we want to emphasize that 
resonant activation is a general phenomenon for phenotypic transitions, which 
are analogous to thermally activated barrier crossing processes. One may find 
application of the idea discussed here to other problems. For example, Spencer 
et al. have shown that cancer cells have persistence behavior similar to bacteria 
|36j . Radiotherapy is a standard cancer treatment option. It normally consists 
of multi-session low dose of radiation in a few weeks. The radiation induces cell 
DNA damage, which eventually leads to apoptosis. In this case it is even easier 
to apply RA. Here radiation is the oscillating signal, apoptosis plays the role of 
antibiotics, and population restoration during the treatment intervals is not a 
serious problem. An optimal strategy may exist on performing the treatment 
utilizing resonant activation. Similar argument applies to chemotherapy. In this 
case stochastic resonance does not exist since the system is not bistable. 
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Figure. 1 A schematic toggle switch model controlling single cell phcno- 
type switch, (a) The toggle switch model. Periodic signal is added to the 
system, perturbing the degradation rate of protein A. (b)Gilliespie simula- 
tions show that the model behaves like a two-state system, corresponding 
to two phenotypes of bacterium. 

Figure 2. Results of single cell simulations, (a) P{tp2n) in the absence of 
the perturbing signal, (b) The mean first passage time (MFPT) versus 
frequency uj of the sine-formed perturbation signal, (c) A weak perturbing 
signal with off-resonance frequency still gives an exponential distribution 
of P{tp2n), here ui = 10hr~^. (d) A weak perturbing signal with resonance 
frequency (red curve) changes P{tp2n) into several separated spikes. Each 
spike overlaps with a peak of the periodic signal (red curve) . uic = 1.6hr~^ . 

Figure 3. Comparison of various strategies at the population level, (a) 
90% quantilc (Q) of Tkui under six different killing strategies, (b) Corre- 
sponding relative antibiotic consumption. Please refer to Table 1 for the 
illustration of killing strategies K. 

Figure 4. Sample population dynamics corresponding to the results in 
Figure 3. (a) Kl. (b) K2. (c) K3, Wc = 1-6 hr-^, T^^ = 27t/u!c = 3.9 hr. 
(d) K4. Black curve shows the dynamics of normally growing population. 
Red curve shows the dynamics of pcrsister population. Gray time win- 
dows in the background indicate antibiotic treatment, blank time windows 
represent environment good for growth. 

Figure 5. Sample population dynamics corresponding to a phenotype 
with slower switching rate (see the main text). Antibiotic treatment is 
either continuous or periodic, (a) Kl. (b) K2. (c) K3, Uc = 0.27 hr~^, 
T^^ = I'KjijJc = 23.3 hr. (d) K4. Black curve shows the dynamics of 
normally growing population. Red curve shows the dynamics of persis- 
ter population. Gray time windows in the background indicate antibiotic 
treatment, blank time windows represent environment good for growth. 



Table 1: Parameters for single cell level simulation 


parameter 


value (reduced unit *) 


notation 




0.05 


basal synthesis rate or A 


7,0 


i 


basal synthesis rate; of B 


^3 




basal degradation rate of A 


^4 


0.56 


basal degradation rate of B 


^5 


5 


dimmer A2 association rate 


k-5 


5 


dimmer A2 disassociation rate 


ke 


5 


dimmer B2 association rate 


k-e 


5 


dimmer B2 dissociation rate 


kr 


5 


binding rate between operon and A2 


k-7 


1 


unbinding rate of A2 from operon 


ks 


5 


binding rate between operon and B2 


k-8 


1 


unbinding rate of B2 from operon 


kg 


1 


synthesis rate of A 


kio 


1 


synthesis rate of B 



*: unit = 130 hr''^ 





Table 2: 


Parameters for population level simulation 


parameter 


value {hr 


notation 


9n 


0.2 


net growth rate of normally growing cells 






under growth condition 


dn 


4 


net death rate of normally growing cells 






under antibiotic condition 


9p 


0.02 


net growth rate of persister cells under 






growth condition 


dp 


1 X 10" 


^ net death rate of persister cells under 






antibiotic condition 



Table 3: 90% quantile of Tkui {Q) under different killing strategies 









Q {hrs) 


notation 








84.6 


no signal + periodic antibiotics 


K2 




^0 


76.0 


no signal + continuous antibiotics 


K3 






33.2 


resonance signal + periodic antibiotics 


Ki 




Eq 


33.0 


resonance signal + continuous antibiotics 


K5 


Sio 


HiSw) 


68.8 


ofF-resonance signal + off-resonance periodic antibiotics 








63.7 


off-resonance signal + continuous antibiotics 



*: H{x) is a Heaviside function which returns 1 if x > 0, otherwise returns 0. 

We assume H = 1 denotes the antibiotic treatment is switched on, and H = 
denotes the antibiotic treatment is switched off. S^^ = Asin{ujct) is the 



sine- formed perturbing signal under resonance frequency Wc- **: Eq represents 
continuous antibiotic treatment (the antibiotic treatment is always kept on). 




Figure 6: P(tp2n) under different signal frequency. All signals are sine-formed, as 
described in the main text. a. oj = 0.3 hr~^ . b. a; = 1.0 hr~^ . c. uj — 2.0 hr~^ . All 
these three cases give values of the MFPT close to the minimum in Figure 2 (b) in the 
main text. 

6 Supplementary Information 

6.1 Large region of u guarantees RA in single cell 
phenotypic transition 

Under resonance frequency Wc « 1.6 hr^^, the mean first passage time 
(MFPT) reaches a minimum. ActuaUy, the MFPT-o; curve displays a 
rather flat bottom, indicating a large region of of the signal frequency 
may generate RA in our system. Figure 6 shows some FPT distribution for 
the p2n transitions (also denoted by P{tp2n)) within this region of signal 
frequency. The P{tp2n) under lo = 0.3, 1.0, 2.0 /ir~^ all produce separated- 
spikes-like distribution. This allows some freedom in real therapeutical 
applications. We have repeated the double well system studied by Schmitt 
et al. [20], and found similar behaviors. 

6.2 Choices of the perturbation target for RA 

In the main text, we reported results with the degradation rate of protein 
A being periodically perturbed. Figure 7 (a-c) shows the MFPT of p2n 
transitions when A's synthesis rate, B's synthesis or degradation rates 
are perturbed at different signal frequency, respectively. Clearly, resonant 




m(hr') 

Figure 7: RA with perturbation on different reactions, a) IVIFPT of p2n transition when 
A's synthesis rate kg is perturbed, b) IVIFPT of p2n transition when B's synthesis rate 
fcio is perturbed, c) MFPT of p2n transition when B's degradation rate is perturbed. 
These perturbations are implemented using sine signal with the same perturbing strength. 
Together with the case of perturbing A's degradation, all four cases generate local min- 
imums around Wc — 1.6hr~^ , but with different fold of change in the MFPT as signal 
frequency varies, d) Stationary distribution on the ua — ns phase plane sampled with a 
long trajectory. Refer to the colorbar for the relative probability of each (nyi,w_g), with a 
decreasing value from the top (red) to the bottom (blue). 
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Figure 8: tp2n under the step-function signal. Other parameters are the same as in 
Table 1 of the main text. 

activation exists in all cases. However, the fold of change between the 
MFPT at the resonant frequency and that with no perturbation varies. 
The effect with perturbation on A's synthesis is similar to that on A's 
degradation, but the effect with either B's synthesis or degradation is 
much weaker. Figure 7 (d) shows the system's stationary distribution on 
the riji — ub plane, which can be related to a potential [37J. Clearly the 
system dynamic shows transition state-like behavior, with the transition 
state having small numbers of both A and B. That is, leaving from the 
state with high A and low B (persister), for the system to make a state 
transition it is more important to reduce the number of A than increasing 
the number of B. This explains why it is more sensitive to perturb the 
reactions involving A. 

6.3 Resonant activation exists with various signal forms 

In the main text, we applied a sine- formed signal to perturb the degra- 
dation rate of protein A in each generic toggle switch, and observed RA. 
A more practical signal form may be unidirectional on the perturbation. 
We tested with a step-function form. That is, the protein degradation 
rate ^3 was only increased from its basal rate periodically by a constant 
value. Again we observed that the MPFT shows a minimum at the same 
Wc as with the sine function form. However, in this case some of the neigh- 
boring spikes in the FPT distribution are not fully separated. Figure 8. 
We observed ~ 9.4% p2n transitions taking place under the growth envi- 
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Figure 9: Distribution of Tkui. a. Case Kl. b. Case K3. Blue dash line shows the 
location of the mean value; red dash line shows the location of 90% quantile of the total 
samples (Q). 



ronment, comparing to ^ 4.3% for the sine function form. Consequently, 
for a system with easy population restoration (i.e., the second example 
discussed in the main text), the sterilization time using the periodic an- 
tibiotics treatment (K3) with the step function signal is longer than that 
with the sine function signal. In this case, strategy K4 under the step 
function signal gives a 90% quantile of the sterilization time ~ 39% of 
that with strategy K2. 

6.4 The distribution of Tkui and the reason for using 
statistics Q 

Figure 9 gives the killing time distribution of the 1000 samples under 
strategies Kl and K3. Clearly both distributions have long tails. Therefore 
we suggest that the 90% quantile may be a better quantity than the mean 
value to compare different strategies, although the conclusions are the 
same. 




Figure 10: Sample population dynamics with corresponding p2n transitions under 
different killing strategies, (a) Kl. (b) K3, ujc = 1.6 /ir"\ T^^ = 2tv/luc = 3.9 hr. (c) 
p2n under Kl. (d) p2n under K3. Black curve shows the dynamics of normally growing 
population. Red curve shows the dynamics of persister population. Gray time windows in 
the background indicate antibiotic treatment, blank time windows represent environment 
good for growth. 




Figure 11: 90% quantile of sterilization time (Q) under signal with different strength. 
Here we use sine-formed signal S = 'ysin{u)t) with strength 7. In the main text, 7 = 
0.15 unit. The maximal 7 should be the basal degradation rate of the target protein, 
which in our case is 0.45 unit. 



6.5 Most p2n transitions take place within antibiotic 
periods under killing strategy K3 

By recording the population dynamics with corresponding p2n transitions 

under different killing strategies, Figure 10 shows that most p2n transi- 
tions take place within the period of antibiotics treatment under strategy 
K3. Therefore, the bacterial population is difficult to re-establish when 
K3 is applied. On the other hand, under Kl, the distribution of p2n is 
exponential, indicating the phenotypic transition may take place during 
either growth condition or antibiotic condition. 



6.6 Signal strength affects the sterilization time of 
bacteria population 

As can be seen from Figure Sll, signals with stronger strength can re- 
duce more time needed for sterilizing bacteria population. However, a 
stronger signal implies a bigger perturbation to the system, which may be 
unfavorable due to practical considerations and toxicity. 
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Figure 12: 90% quantile of sterilization time (Q) under different duration of antibiotic 
treatment. 

6.7 Duration of each session of antibiotics treatment 
affects sterilization time 

At resonance frequency, the p2n transition distribution is multi-spike shaped, 
This allows discrete antibiotics treatment, as long as each session covers all 
or most of the p2n transition time. We performed population level simu- 
lations with different antibiotics treatment duration centered at the signal 
peak times. Figure 12 shows the 90% quantile of the sterilization time. As 
expected, upon increasing the antibiotics treatment duration from zero, 
the sterilization time first drops sharply, and then fiattens. Further in- 
creasement of the antibiotic duration does not decrease the sterilization 
time. Therefore in this case one may further reduce the duration of each 
session of the antibiotics treatment to ^ 1.3 hr from the 2hr value used 
in in the main text, without significant increase of the sterilization time. 
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